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ABSTRACT 

We consider constraints on the structure formation model based on non-Gaussian 
fluctuations generated during inflation, which have \ m distributions. Using three data 
sets, the abundance of the clusters at z = 0, moderate z and the correlation length, 
we show that constraints on the non-Gaussianity and the amplitude of fluctuations 
and the density parameter can be obtained. We obtain an upper bound for Q m and 
a lower bound for the non-Gaussianity and the amplitude of the fluctuations. Using 
the abundance of clusters at z ~ 0.6, for the spectrum parameterized by cold dark 
matter (CDM) shape parameter Y = 0.23, we obtain an upper bound for the density 
parameter Sl m ~ 0.5 and lower bounds for the amplitude a% ~ 0.7 and for the non- 
Gaussianity of fluctuations G ~ 2 (m <~ 200), where G = 1 for Gaussian. 
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1 INTRODUCTION 

It is generally accepted that the primeval density fluctua- 
tions have grown gravitationally to form the present struc- 
ture of the Universe. Many models based on Gaussian adia- 
batic fluctuations generated during inflation have been dis- 
cussed. However, the standard CDM model, i.e. the struc- 
ture formation model in a spatially flat universe filled with 
cold dark matter which has Gaussian adiabatic fluctuations, 
suffers from several problems. One of the most severe prob- 
lems is the lack of the cluster abundance at moderate z. This 
is caused by the fact that in a flat universe filled with the 
matter, fluctuations grow rapidly. To overcome this problem, 
a low-density universe model, with or without cosmological 
constant has been investigated extensively. 

Recently an alternative was proposed, a model based on 
non-Gaussian isocurvature fluctuations produced by mas- 
sive scalar fields frozen during inflation, which have Xm dis- 
tributions (Peebles 1997a,1999a,b). The most notable point 
of this model is the early formation of objects because of 
its non-Gaussianity. If the probability distribution function 
(PDF) has a broad tail compared with a Gaussian PDF, 
the evolution of clusters becomes slow. In this model there 
are three important parameters; the amplitude and the non- 
Gaussianity of the fluctuations and the background density 
parameter. It is important to explore methods of constrain- 
ing these parameters. 

It is recognized that the non-Gaussianity affects the cor- 
relation length of clusters. Robinson, Gawiser & Silk (1998) 
suggest that by combining the observations of the abundance 
of clusters at z = and the correlation length, constraints on 



the amplitude and the non-Gaussianity of the fluctuations 
are obtained simultaneously for a fixed density parameter 
O m . As mentioned above, because the density parameter 
and the non-gaussanity affect the evolution of abundance, 
the observations of the abundance of clusters at moderate z 
give another strong constraint on the model. Much work on 
the effect of changing Q m on the evolution have been carried 
out assuming Gaussian initial conditions (Frenk et al. 1990; 
Eke et al. 1996; 1998; Viana & Liddle 1996,1998; Fan et al. 
1997; Bahcall & Fan 1998; Henry 1997; Carlberg et al 1997). 
Little work has been carried out, however, on the effect of 
the non-Gaussianity, except for a topological defect model 
(Chui & Ostricker 1998). 

In this paper, we extend Press-Schechter (PS) theory to 
non-Gaussian initial conditions and make clear how the non- 
Gaussianity affects the evolution of abundance and the bias. 
Then we show that combining the observations of the abun- 
dance of clusters at moderate z with the method proposed 
by Robinson. et.al, we can obtain a constraint on background 
cosmology and fluctuations simultaneously: an upper bound 
for Q m and a lower bound for the non-Gaussianity and the 
amplitude of the fluctuations. 



2 MASS FUNCTION AND BIAS: FOR 

NON-GAUSSIAN INITIAL CONDITIONS 

We extend PS theory (Press & Schechter 1974; Bond et al. 
1991) to non-Gaussian initial conditions and calculate the 
mass function and the bias of clusters. 

In PS theory, the region of the scale R with overden- 
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sity greater than a critical amount S c is collapsed to form 
a bound object of mass M = (47r/3)p-R 3 , where p is the 
mean density of the universe. The threshold is assumed to 
be given by the spherical collapse model. For Gaussian initial 
conditions, it has been shown to be a good approximation 
(Lacey & Cole 1993). We assume that this approximation is 
good for non-Gaussian initial conditions. The probability of 
a region with mass M collapsing at redshift Z\ is given by 

f°° 8 
F(8 1 ,M) = / P{S)dS, 5 1 = T ^- T , (1) 



D(zx)> 



where P(S) is the PDF in question and D(z) is a linear 
growth factor. According to Jedamzick, the above quantity 
can be written as 

1 f°° 

F(5i,M) = ~ Mn(M',zi)P(M\M)dM', (2) 
P hi 

where n(M, Zi) is the mass function at z = Zi and P(M\M') 
denotes the probability of finding a region with mass M 
overdense by Si or more, provided it is included in an iso- 
lated overdense region with mass M' (Jedamzick 1995). We 
assume this quantity is given by 



P(M\M') = P(8 M > 8i\5 M < = h) 



P{S\6i)dS=l/f,{3) 



where P(5\8i) is a conditional probability function. Hence, 
the mass function n(M, zi) is given by 

fp dF(S u M) 



n(M,zi) = 



(4) 



M dM 

The bias of the bound objects can be calculated using 
the formalism developed by Mo and White (Mo & White 
1996). Consider the collapsed region of scale Ri with over- 
density 8i and the uncollapsed region of Ro with So- If at 
present <5i = 8 C /D{z\), then the halo is identified as being 
formed at Z\. The number of objects of Mi identified at Z\ 
in a spherical region with a comoving radius Ro is 

,M d 



JV(l|0)dAfi = -/ 



Mi dMi 



dSP(8\8o) 



(5) 



<5i 



Thus, the average overdensity of objects in the sphere of ra- 
dius Ro relative to the global mean halo abundance becomes 

JV(1|0) , 



%(1|0) = 



(6) 



n(M, z)Vo 

in Lagrangian space (Mo & White 1996). 

If the region is sufficiently large compared with region 
1 ( Ro ~> Ri, {81} > {8q}), regardless of the statistics, the 
relation 



<W- p > = (SS) 



(7) 



is satisfied (Appendix A) and the conditional probability is 
given by 

P(8i\8o) =P(8i -So). (8) 

(Appendix B). If this is the case, we obtain the normaliza- 
tion of the mass function 



/ = !/(/ P(S)d8), 



(9) 



which agrees with the one introduced to ensure that all the 
mass in the universe is accounted for (Lucchin and Matarrese 
1989, Chui and Ostriker 1998). Since 



JV(l|0)dM = -/ 



Mo d 
Mi dMi 



dSP{8), 



(10) 



the overdensity of the bound objects is given by 



ra(l|0) - n{M) 
n(M) ' 



n(l|0)=n(M)| 4l ^, 1 _, . (11) 



This reproduces the peak-background split formalism 
(Kaiser 1984; Bardeen et al. 1986; Robinson et al. 1998). 
Furthermore, we assume that reduced moments 



Ajv = 



(s N ) 

{8 2 ) N / 2 



(12) 



do not depend on mass M on the relevant scales. 

In the following discussion, we make these assumptions 
and use above formalisms. Note that the accuracy required 
to use Press-Schechter theory for non-Gaussian initial con- 
ditions is currently uncertain. One way to confirm its valid- 
ity is to perform N-body simulations. Recently Baker and 
Robinson have carried out N-body simulations for several 
non-Gaussian models to show mass function of PS formal- 
ism agrees well with that observed in simulations (Robinson 
& Bakerl999). It will be necessary to perform further de- 
tailed tests of the accuracy of this formalisms, especially for 
the bias model. 



3 THE ABUNDANCE AND THE BIAS OF THE 
BOUND OBJECTS. 

We use the formalisms developed in the previous section and 
investigate the effect of non-Gaussianity on the evolution 
of abundance and the bias of the bound objects. First we 
consider the evolution of the number density. We define the 
rate of change with z of the number density by 



An(M, z) = 



dn(M, z) 



n(M, z) dz 
This can be written as 



An(M,z) = 

where 
v{M,z) 



dlog vP(y) 



dv 



u{M,z) 



dlogD(z) 
dz 



(13) 



KM, z),(14) 



<t(M)D{z) 



*(M)=. 8 (-) 



M \ ~( 3 + n )/ 6 



(15) 



where M* is the mass contained in the spherical region of 
8ft -1 Mpc and n is the spectrum index. This indicates that 
for a fixed PDF, if fluctuations decrease slow with z (as in 
low Q m universe), the number density also decreases slowly 
with z. We fix the background cosmology and compare the 
evolution of the number density for different PDFs. The dif- 
ference of this quantity between a Gaussian PDF and the 
PDF in question can be written as 



Ana — An — 



d\ogg(u) 



dv 



v{M) 



d\ogD(z) 
dz 



u{M). (16) 



Here g(v) denotes the ratio of the PDF to a Gaussian PDF 
P(y) 



(17) 
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where Pg is a Gaussian PDF. Next we calculate the bias 
of clusters at the present time. The linear bias 8 c i = bS in 
Eulerian space is given by 

dlogn(M) 



b(M) 



d5 c 



The difference from a Gaussinan PDF is given by 



b G (M)-b(M) = 



d\ogg{v) 



dv 



where the bias for a Gaussian PDF is 
b G (M) = 1 



o{M) ' 



v{Mf - 1 



(18) 



(19) 



(20) 



(Cole & Kaiser 1989; Mo & White 1996). 

If a PDF has a broad tail compared with a Gaussian 
PDF, therefore, the number density of the objects M evolves 
slowly with z and the bias of such objects is small compared 
with those for a Gussian PDF. This is natural because these 
rare objects can be formed easily for a PDF with a broad 
tail. 



4 THE x 2 MODEL 

In this section, we investigate the Xm model ( Coles and Bar- 
row 1987, Moscardini et.al 1991, Weinberg and Coles 1992). 
In this model, fluctuations are drawn from Xm distributions. 
Peebles gives one realization of this model in the framework 
of inflationary universe(Peebles 1997,1999a). In his model, 
CDM fields are squeezed massive scalar fields which have 
O(N) symmetry. The isocurvature perturbations are gener- 
ated during inflation. The density and the over-density of 
CDM fields are given by 



Pcdm 



(21) 



where m is the number of the CDM fields <j>, which have 
0(m) symmetry and a Gaussian PDF and fi is the mass of 
the CDM fields. Then, the PDF of the overdensity is given 
by 



(1 + /T^m/a-l 




dv.{22) 



If the autocorrelation function of the Gaussian field (p(x ) is 
£^ = (4>(x 1 )(j)(x2)} oc Xi2(xi2 = \&i — x 2 \), the autocor- 
relation function and the reduced moments of 5 are given 
by 



(*(xi)*(x 2 )) 



2 > 



Aw = 



2\N/2 1 (( Xl2X23 . ■x N - 1N XNl) C ) 



5 - a «>W2) 



(23) 



where angular brackets denote averaging on some smoothing 
length (Peebles 1999b). 

Note that two problems seem to exist in this model. 
First, because the autocorrelation function is positive for all 
separations, the integral of £(r) over all space does not van- 
ish (White 1998). Next it can be shown that the reduced mo- 
ments remain approximately same for a variety of smoothing 



scales if e is scale-invariant (Peebles 1999b, White 1998). This 
seems to be in contradiction to the central limit theorem 
that states all distributions must become Gaussian when 
smoothed on sufficiently large scales. If desired, however, it 
is possible to modify the Peebles model so that these diffi- 
culties disappear. For example, an off-centered Xm model, 
where p — /x/2 — 4>o) 2 , produces a distribution that is 

non-Gaussian on small scales and Gaussian on large scales. 
This scale is determined by the mass scale <f>o- Here, we are 
interested in the scale relevant to cluster formation where £ 
has a power-low spectrum for which £ > and the smooth- 
ing length is not so large, so we believe that these technical 
difficulties do not affect us. Indeed, our method can be ap- 
plied to any non-Gaussian model on this scale. 

In the Xm model the non-Gaussianiy is represented by 
one parameter m. As the effect of the non-Gaussianity de- 
pends on the broadness of the PDF, we introduce another 
parameter G to represent the amount of deviation from 
Gaussianity: 



G = V2^ I P{y)dv/ I e 



dv, 



(24) 



which is the likelyhood relative to Gaussian of 3a or rare 
events (Robinson et al. 1998). G is related to m as follows: 
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5 EVOLUTION OF ABUNDANCE AND 

CORRELATION LENGTH OF CLUSTERS 

We examine the evolution of the number density of rich 
clusters from z = to z = 1. We consider the clusters 
with the mass M(< R com = 1.5/i _1 Mpc) > 8 x 10 14 M Q , 
where M(< R CO m) represents the mass contained within the 
comoving length R CO m (Bahcall & Fan 1998). The number 
density of such clusters is given by 

N(M(< R com = 1.5/i _1 M P c) > 8 x 1O 14 M ) = / n(M)dM. 

J M lim 

We relate the mass threshold M(< R com = 1.5h~ Mpc) = 
8 x 10 14 Mq to the virial mass Mu m using the density profile 
M(< R) oc 7? 64 (Bahcall & Fan 1998). The result is shown 
for different Q. m (G = 1, A = 1 — Q m ) and G (Q m = 1) 
in Fig.l. The Gaussian model in the universe of Q m = 1 
is strongly disfavored (Calberg et al 1997; Bahcall & Fan 
1998). The non-Gaussianity of fluctuations makes the evo- 
lution of the abundance slow. 

The correlation length ro of clusters with the mean sep- 
aration d = 7V(> i?)~ 1//3 is obtained as follows (Robinson et 
al. 1998). Using the mass function we obtain the scale R of 
clusters with the mean separation d = N(> i?) _1//3 . Then, 
the amplitude of two-point correlation function of matter 
fluctuations at correlation length is given by £ m = l/b(R) 2 . 
If we know and the spectrum index, n, we can obtain the 
correlation length. The bias is given by 



1 + 



1 



L V V m ' 



(25) 
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Figure 1. The evolution of cluster abundance of M(< R CO m = 
1.5/i _1 Mpc) > 8 X 1O 14 M clusters for a different Q m (G = 1, 
A = 1 — Q m , upper panel) and for a different G (Q m = 1, lower 
panel). 

For large G (small m), the bias becomes small and the corre- 
lation length becomes short. We fix erg = 0.5, n = — 1.4 and 
Q m = 1; the bias and correlation length of the clusters with 
mean separation d = 30ft -1 Mpc are then given as follows: 
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4.02 


2.04 
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6 CONSTRAINTS ON Q M , G AND cr 8 

We can give constraints on Q m , G and as by using the abun- 
dance of the rich clusters at z — 0, z = 0.575 and the corre- 
lation length. We use the following data sets. 

• The abundance of M(< R com = 1.5) > 8 x lO 14 ft _1 M 
clusters at z=0 is taken from temperature function of T > 
6.4 kev (Bahcall & Fan 1998). We use 2±\ X W^Mpc" 3 
(Henry & Arnaud 1991). 

• The abundance at z=0.575 is taken from the data of 
the Einstein Medium Sensitivity Survey (EMSS). We use 



n(z = 0.5 - 0.65) = 3±?;g x W~ 8 h 3 Mpc -3 for Q. m = 1. 
The error bars represent 2a statistical errors (Bahcall & Fan 
1998). (For fi m < 1 we use n(z = 0.5-0.65) = 2t\ x 10~ 8 ft 3 
Mpc" 3 .) 

• The correlation length is taken from the Auto- 
matic Plate Measuring (APM) Survey. We use ro = 
14.2l^/i _1 Mpc for the clusters of mean separation d = 
30/i _1 Mpc. The error bars represent 2a statistical errors 
(Croft et al. 1997). 

Henceforth, for simplicity, we only investigate the model 
with a cosmological constant and use the spectrum param- 
eterized by CDM shape parameter F (Bardeen et al. 1986). 
We will assume F = 0.23 (Vianna & Liddle 1998). 

From the above data sets, we can give constraints on 
Qm, G and a$ (Fig. 2). First, consider the abundance at 
z = 0. For large G, the clusters are formed easily so we need 
small as- Since (Tg e // = ag£lm +3 ^ 6 determines the abun- 
dance of clusters, we need large as for small O m . The fit 
contour runs from the upper left (large G and small erg) to 
the lower right (small G and large as) in the as vs G plane 
and this moves to right as Q m becomes small. Next, consider 
the abundance at z ~ 0.6. For fi m ~ 1, if fluctuations are 
Gaussian (G = 1), the abundance at z ~ 0.6 is too small 
compared with the data of the EMSS. For large G, because 
of the effect of the non-Gaussianity, the abundance becomes 
consistent with the observation (Fig.l). The fit contour runs 
in the plane in the same direction as of z — 0, but the slope 
is smaller. As fl m decreases, the abundance at z ~ 0.6 can 
be explained even if fluctuations are Gaussian. Then, as fl m 
decreases, the fit contour moves to right more slowly than 
that of z — 0. The region that can support these two data 
sets at the same time therefore exists in the upper left of the 
plane for large Q. m and stretches to the lower right as Q m 
becomes small. 

Next consider the correlation length. At the scale of 
the correlation length, the variance of matter fluctuations is 
given by a = 1/b. For large G, therefore, because the bias 
6 is small, we need large as- The fit contour runs from the 
lower left to the upper right of the plane, i.e. orthogonal 
to the fit contours of the abundance (Robinson et al. 1998). 
The correlation length is mainly determined by the statistics 
of fluctuations, which does not change as Q m changes. The 
region that supports the correlation length and abundance 
at z — therefore moves from lower left to upper right as 
Q m decreases. Now, we combine all the data. For O m ~ 1 
we cannot find a region that will support all the data sets. 
As O m decreases, we come to some Q m , for which the region 
that supports all the data sets exists. This region moves to 
the upper right as fi m becomes small; we can then obtain 
the upper bound for fi m and the lower bound for os and G. 

From the above data sets, the upper bound Q m < 0.5 
and the lower bound G > 2 and as > 0.7 are obtained. 
Particularly in the universe of Q m = 0.3, which is favored 
by other observations, non-Gaussianity of the order G ~ 4 
is favored. 

The shape of the cluster temperature function could 
possibly give another constraint on the model (Kitayama 
and Suto 1996). We show the temperature function at the 
present time for the parameters that are consistent with all 
three data sets (fi m = 0.3, A = 0.7, G = 4.02, as = 0.75 and 
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Figure 2. Fit contours of the three data sets in trg vs G plane. 




10° 10 1 
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Figure 3. The cumulative cluster temperature function for the 
X^j model at the present time. 



Q m = 0.1, A = 0.9, G = 6.33, a 8 = 1.1 (Fig.3)). Here we 
convert the virial mass M v to X-ray temperature using 

where A„ is the virial over-density (A„ = Pviriai/pc) and 
we take f3 = 1 (Eke et al. 1996; 1998; Borgani et al 1998). 
The shape of the temperature function agrees well with the 
data of Henry and Arnaud for high temperatures (Henry & 
Arnaud 1991, Henry 1997). Note that as Q m decreases, the 
non-Gaussianity and ag becomes large, so the slope of the 
temperature function becomes smaller. Then, the difference 
in the abundance for different fl m is larger for low temper- 
atures. This provides us with the possibility of obtaining 
another constraint on the low-density universe model. 



7 CONCLUSION 

We have given constraints on the Xm model, where fluctu- 
ations are drawn from Xm distributions. We extended the 
Press-Scheter theory to calculate the abundance and bias of 
clusters. The non-Gaussianity of fluctuations makes the evo- 
lution of the abundance slow. The model with non-Gaussian 
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fluctuations of the order G ~ 16 (m = 1) in Q m = 1 universe 
has roughly the same evolution of the number density as a 
Gaussian model in a Q m ~ 0.2 universe, where G represents 
the non-Gaussianity of the fluctuations (G = 1 for Gaus- 
sian) . On the other hand, the strong non-Gaussianity of fluc- 
tuations make the correlation length too short. Combining 
the three data sets, the abundance of the clusters at z = 
and z ~ 0.6 and the correlation length, constraints on the 
non-Gaussianity and the amplitude of fluctuations and the 
density parameter have been obtained. We have shown that 
the upper bound for fi m and the lower bound for the non- 
Gaussianity and the amplitude of the fluctuations can be 
given. For the spectrum parameterized by CDM shape pa- 
rameter F = 0.23, we have obtained an upper bound for the 
density parameter fi m ~ 0.5 and lower bounds for the am- 
plitude as ~ 0.7 and the non-Gaussianity G ~ 2 (m ~ 200) 
of fluctuations. In the universe with Q. m = 0.3 and A = 0.7, 
non-Gaussianity of the order G ~ 4 (m ~ 30) is preferred. 
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APPENDIX A: CROSS-CORRELATION 
BETWEEN SCALES 

In this appendix, we show that if Ro^$> Ri, the relation 
(W P > = W>. (Al) 

is satisfied. This is the extension of the relation for Gaussian 
initial conditions which was derived by Bower under the 
assumption that the volume of the region is sufficiently 
larger than that of the region 1 (Bower 1991). The smoothed 
overdensity in the region of scale R is 

S(x,R) = J 0^5 k W R (kR)e-* k x , (A2) 

where Wn(kR) is a smoothing function. The n-point func- 
tions of the random field 5^ are expressed as 

( S k 1 S k 2 ) = £a(fci,fca), 
(hjkjkj = 6(fci,fe,fc 3 ) 

(Sk^^kJkJ = £4(fcl J fc2 1 fc3,fc4)+&(fcl,fc2)&(fc 3 ,fc4) 

+cycl. (2-terms). (A3) 

Consider the region 1 with scale Ri centered xi contained 
in the region with scale Rq centered xq. The variance of 
the overdensity in the region i is 

a 2 (i?0 = (5?) = J J^LwUkRi)P(k), (A4) 

where P(k) is defined by £ (fej, kj) = 5o(ki + kj)P(ki). The 
covariance is written as 

a 2 (R R 1 ) = (8{R )5{R 1 )) 

-^ s P(k)W(kR )W(kR 1 )(e- i fe ' y )i e0 , (A5) 

where y = Xo — xi. Averaging over y is done as follows: 
select a random set of points y for which all points in the 
region 1 are independent and enclosed in the region 0, then 
average over this random set of y. The averaged term can 
be written to appear like another window function 



W R (kR) 



)ieo 

d 3 y w(y)e~' k y / I d 3 yw{y), 



(A6) 



where R is the scale covered by the points y for which re- 
gion 1 lies in region 0. The weight function w(y) is defined 
so that every point in the region has an equal chance of 
contributing to the average density measured in the region 
1. Then we can write the real-space window function of the 
region as the convolution 



W Ro (x) 



d 3 x'W Rl (x)W R (x - x) 



(A7) 



In Fourier space 
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W Ro (kR ) = W k {kR)W Rl {kRi 



(A8) 



Note that it is possible to find real weighting function only 
if the region 1 is sufficiently small compared to the region 0, 
i.e. Ro 3> Ri- If this is the case, the covariance is given by 



a(R Ri) 



j£k_ 

(2tt) 3 



WUkRi)P{k) = g\R ). 



(A9) 



Extensions to the higher moments are quite similar. For ex- 
ample, 



d 3 fci ■ ■d 3 fe n <5| ) (fci + ■ ■ kn)(6 ki ■••S kn ) 



xW Ro (kiR ) ■ ■ ■ W Ro (k n - p R ) 

X [W Rl (kn-p+lR^W^kn-p+lR)} 

x[W Rl (k n Ri)W n (k n R)] 

' 3 U J 3 1 c3 



(Tfci ■ -(Tk n 5h{ki + ■ ■ k n ){S ki ■ ■ ■ S k J 

xWRgi^Ro) ■ ■ ■ W Rg (k n R ) 
= (So)- (A10) 

APPENDIX B: CONDITIONAL PROBABILITY 

We will derive the Edgeworth expansion (Bernardeau & Kof- 
man 1995; Juszkiewicz et al. 1995) for the conditional prob- 
ability of So and Si in the case that Ro 2> Ri is satisfied. 
The generating function for joint moments 



C(/z, u) = \ dSdGe 



ifj,5-\-iiy6 



P(5,6) 



gives rise to 

d n+m C(n,u) 



(s n e m ) = 



d(iy) n d(iv) r 
The characteristic function becomes 



K{n, v) = log C{n,v) 



E 

(n,m)#(0,0) 



~ Anri 



where connected joint moments 

A mn = (8 n 6 m ) c 



(Bl) 



(B2) 



(B3) 



(B4) 



are used. The inverse transformation gives the conditional 
probability function 



p(6,e) = 



i 



(277) 



dfidue 



-i/j,5 — iv6 



(2?r) 2 



dfidue 



-iflS—iu0— 77 /i. 2 — 7f f 2 —rfiu Yi(fi,v) 



where 

E(/i, !/) 



E E 

JV=3 m + n=]V 



(B5) 



(B6) 



is the nonlinear part of the characteristic function and r — 
(59) is a cross-correlation. Here we rescale S and 6 so as to 
have a unit variance. Simple calculation yields 



P(S,0) = e E(i ^' 1 



1 (s- r e)^ 1 2 



27r\/r 



(B7) 
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To proceed further, we use the relation derived in appendix 
A: 



(S N ) = Ajvo, {S p e N ~ p ) = r p \ N, (9 N ) = Xon. 
Now, let us define 
5 -rG 



vT 



We then obtain 

/ 00 

P(8,6) = exp [J2 



(B8) 
(B9) 



VI - r 2 da 



1 _I Q 2 

x^=e 2 



2tt 



where 
P(0) = exp 



E 

v iV=3 



Aqjv 



1 

"M 



(BIO) 



2tt 



2tt 



e 2 



l + ^3(9)+0(AU 4 ) 



Here we use the definition of the Hermite functions 
d 



e 2 



Then the conditional probability function is given by 
P(S,6) 



(Bll) 
(B12) 



P(8\9) = 



P(9) 



exp 



E 

\N=3 
1 -la 2 



■^■(Ajvo — Aojvr^) 



8 



VI - r 2 9a 



2tt 



(B13) 



Consider the small region 1 contained in a large region 0. 
Defining the variance o\ for <5i and oo for So and rescaling 
yields the Edgeworth formula for conditional probability 



P{Sx\So) = 



27ryV| - a 2 



exp 



(Sr-So) 2 



J, (^)-(^o 3 ) g f ft-&> V 

Using the relation 

(5?) >W, 

the conditional probability is reduced to 
P(*i|tfo) = - So). 



(B14) 

(B15) 
(B16) 



